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Abstract 

Impact mixed-mode I/II fracture experiments were carried- 
out on rigid foam using single edge notched bend (SENB) 
specimens with inclined notch. A drop-weight rig was used 
to perform the low-velocity impact fracture tests. SENB 
specimens with different notch inclination were tested at 
various impact potential energies by varying the drop-height 
in order to determine the critical impact potential energy at 
which the crack growth is initiated from the notch tip. Three- 
dimensional linear-elastic finite element simulations of 
impact behavior of the SENB specimens were carried-out. 
The virtual crack closure technique was applied to analyze 
the strain energy release rate under impact loading. The 
impact fracture toughness was characterized in terms of the 
critical dynamic strain energy release rate at the onset of 
crack growth. It was found that the critical dynamic strain 
energy release rate decreases with increase of notch 
inclination. This finding was attributed to the decrease of the 
mode II component of the loading. SENB specimens with 
different notch inclination were tested also under static 
loading. The static fracture toughness was extracted from the 
test data by the virtual crack closure technique in 
conjunction with finite element simulations. It was found 
that the impact fracture toughness of the rigid foam 
investigated is between 22.4 % and 32 % lower than the static 
one at various notch inclinations. This finding indicates that 
the static fracture toughness should not be used in the 
fracture mechanics based design of rigid foam structures 
which are subjected to impact loading during their lifetime. 
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Introduction 

Technical development facilitates an extensive use of 
foam core sandwich structures for load-bearing 
components in modern engineering (shipbuilding, car 
industry, aeronautics), where high strength and 
stiffness per weight are of primary importance 



(Gibson and Ashby, 1988; Reyes and Rangaraj, 2011; 
Christensen, 2000; Quispitupa et al., 2009; Shi-Dong 
Pan et al. 2008; Rangaraj, 2010). A sandwich structure 
typically consists of two thin stiff face sheets bonded 
to a relatively thick lightweight core. The face sheets 
mainly carry bending moments and axial loads as 
compressive or tensile stresses, while the core works 
in shear supports the faces against bucking and 
wrinkling. Sandwich structures provide additional 
advantages like very well acoustic and thermal 
insulation, high corrosion resistance, competetive 
energy absorption and buoyancy. The ability of a foam 
core sandwich structure to sustain loads essentially 
depends upon the properties of the foam core. 

Crack initiation and propagation in the core 
significantly reduces the stiffness and may lead to 
failure of the sandwich structure. Therefore, it is of 
great practical importance to develop experimental 
and theoretical methods for investigation of the rigid 
foam fracture behavior. 

Static fracture behavior of cellular foams has been 
previously investigated by Ashby et al. (1984). Static 
fracture in closed-cell foam under mode I loading 
conditions has also been studied by Bazant et al. (2003). 
Mode III static fracture toughness of rigid foam 
materials has been characterized by Farshad and 
Flueler (1998) using an anti-clastic plate bending test 
method. A special loading device has been applied for 
inducing almost a pure mode III crack loading 
conditions in the foam specimen. An experimental and 
numerical investigation of debonding in foam core 
sandwich structures under static loading has been 
carried-out by Prassad and Carlsson (1994). Fracture 
tests have been performed using double cantilever 
beam and shear sandwich specimens. Test data have 
been analyzed by linear-elastic fracture mechanics 
concepts. For this purpose, two-dimensional finite 
element models have been developed. The models 
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have been verified by comparisons between the 
calculations and the measurements. Goswami and 
Becker (2001) have analyzed a static delamination 
crack along the face sheet/core interface in a sandwich 
structure subjected to transverse loading by a finite 
element model. The effect of different structural 
parameters on the face sheet/core delamination 
fracture has been studied. The ParisTaw equation has 
been applied to characterize fatigue crack growth in 
foam material by Noble and Lilley (1981). The 
environmental effects on the fatigue fracture behavior 
of Polycarbonate foam material have been analyzed by 
Yau and Mayer (1986). 

The present article describes an investigation of low- 
velocity impact fracture behavior of rigid foam using 
SENB specimens with inclined notch (it should be 
mentioned that such investigation has not been found 
in the literature available). In the experimental part of 
the study, low-velocity impact fracture tests were 
performed on SENB specimens with different notch 
inclinations. The basic purpose of the testing was to 
determine the critical impact potential energy and the 
fracture time at the onset of crack growth for different 
mixed-mode I/II crack loading conditions. Three- 
dimensional linear-elastic finite element simulations of 
the SENB specimens were carried-out in order to 
calculate the critical dynamic strain energy release rate 
using the experimental data. 

Experimental Procedure 

The foam under investigation in the present article is 
Divinycell HCP100. This is polyvinylchloride closed- 
cell foam with nominal density of 400 kg/m 3 . 

Steel cylindrical impact™ 
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FIG. 1 GEOMETRY AND LOADING OF THE SENB SPECIMEN 
(THE ANGLE OF INCLINATION IS MARKED BY a) 

Single edge notched bend specimens were cut-out by a 
diamond blade saw from foam panels according to 
FIG. 1. The length and the height of the specimen were 
L=220 mm and H=50 mm, respectively. The specimen 
width, W, was equal to the foam panel thickness (30 



mm). An inclined notch with length of 25 mm was 
introduced as initial crack in the specimen mid-span. 
The notch inclination was varied from 15° to 90° with a 
step of 15° for the different specimens. The notch tip 
was sharpened by pressing a new razor blade in it 
prior to testing. None of the specimens was pre- 
cracked. 

The SENB specimens were supported on two 12 mm 
diameter steel cylinders. The support span length was 
200 mm (FIG. 1). The impact fracture tests were 
performed using a drop-weight rig. During testing the 
specimens were impacted at the mid-span by a steel 
cylinder with a diameter of 18 mm and length of 40 
mm. After the first impact, the steel cylinder was 
captured to avoid secondary impacts. The impacting 
mass, 4.13 kg, was kept constant for all the tests. A 
range of impact potential energies were attained by 
varying the drop-height at each notch plane 
inclination in order to determine the critical energy. 
The tests were carried-out at room temperature (23°C). 
Five replicate specimens were tested for each notch 
plane inclination. High repeatability of the test data 
was obtained. A computer system was used for data 
acquisition during testing. The contact force-time 
response was monitored continuously at each test. 
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FIG. 2 EXPERIMENTALLY OBTAINED CONTACT FORCE-TIME 

RESPONSES OF THE SENB SPECIMEN WITH NOTCH 
INCLINATION OF 15" AT IMPACT POTENTIAL ENERGY OF 9.6 J 
(CURVE 1), 13 1 (CURVE 2), 14 ] (CURVE 3), AND 14.8 J (CURVE 4). 
THE CRACK GROWTH ONSET IN CURVE 4 IS SIGNALLED BY A 
SUDDEN DROP. THE TIME OF CRACK GROWTH INITIATION 
(I.E., THE FRACTURE TIME) IS tf =5.5 ms 

The onset of crack growth from the notch tip was 
determined by inspection of the specimen lateral 
surface with a microscope after each test (the notch tip 
position was marked on the specimen surface prior to 
testing). In this way it was found that the crack growth 
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in SENB specimen with notch inclination of 15° was 
initiated at impact potential energy level of 14.8 J. It 
should be noted that SENB specimens with notch 
inclination of 15° were tested at different impact 
potential energies. No crack growth was observed at 
energy levels lower than 14.8 J. Contact force-time 
curves for the various impact potential energies used 
in SENB specimen with notch inclination of 15° are 
depicted in FIG. 2. The data in FIG. 2 indicate that the 
maximum contact force generated during an impact 
fracture test increases with the impact potential energy, 
while the duration of the impact event (i.e. the time for 
which the impactor and the specimen are in contact) is 
almost constant for different impact potential energies. 
It was found that the crack growth initiation was 
signaled by a sudden drop in the contact force-time 
curve for the critical impact potential energy, Epc, 
(FIG. 2). This finding was attributed to the degradation 
of the SENB specimen bending stiffness due to the 
crack growth initiation. The time at the drop was 
recorded as the fracture time, tf. It can be summarized 
that the experimental investigation yielded two 
important parameters (critical impact potential energy, 
Epc=14.8 J, and fracture time, tf=5.5 ms) which were 
used in the finite element analysis of the SENB 
specimen with notch inclination of 15° in order to 
evaluate the impact fracture toughness. 

TABLE 1 MEASURED CRITICAL IMPACT POTENTIAL ENERGY, 
EPC, AND FRACTURE TIME, tC, AND CALCULATED AVERAGE 
TOTAL DYNAMIC CRITICAL STRAIN ENERGY RELEASE RATE, 

GDa /~i Da / /~i Da /~* Da / /~i Da 

c , MIXED-MODE RATIOS, IC C IIC C / AND, 
fjDa i Da 

IIIC C , AND AVERAGE TOTAL STATIC CRITICAL 

STRAIN ENERGY RELEASE RATE, c , IN SENB SPECIMENS 
WITH DIFFERENT NOTCH PLANE INCLINATION, a. 
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Epc, J 


14.8 


13.4 


12.2 


11.5 


10.8 


10.4 


tc, ms 


5.5 


5.6 
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6.6 
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C ,I/m 2 


3970 


2853 


2343 


1907 
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176 


s-*Da i /~*Da 

u IC 1 u c 


0.51 


0.53 


0.59 
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0.81 


1.00 


/~i Da / s~i Da 

'-'iic 1 °c 


0.43 


0.41 


0.35 


0.33 


0.17 


0.00 


/-i Da 1 s-tDa 
^IIIC 1 u c 


0.05 


0.04 


0.04 


0.02 


0.00 


0.00 


G S r a 


5117 


3908 


3301 


2801 


2492 


224 



The experimentally measured values of critical impact 
potential energy level and fracture time at notch 
inclination varying from 15° to 90° with a step of 15° 



are summarized in TABLE 1. One can see that the 
critical impact potential energy decreases with 
increasing the notch inclination. This finding was 
attributed to decrease of mode II component of crack 
loading and to decrease of SENB specimen bending 
stiffness with increase of the notch inclination. 

Test Data Analysis 

The impact fracture toughness was obtained from the 
test data using methods of linear-elastic fracture- 
mechanics. For this purpose, three-dimensional linear- 
elastic finite element models of the SENB specimens 
with different notch inclination under impact loading 
were developed using the ABAQUS/Explict program 
system. 




FIG. 3 FINITE ELEMENT MESH USED IN THE SIMULATIONS OF 
THE SENB SPECIMEN WITH INCLINED NOTCH 

The geometry and dimensions of the models were 
dictated by the SENB specimen configurations used in 
the low-velocity impact fracture tests. Due to the fact 
that the specimen had a plane of symmetry, and the 
width of the finite element model was reduced to 15 
mm (i.e., only half of the SENB test specimen was 
considered in the modelling). This was done in order 
to reduce the finite elements number in the model and 
to decrease the computational time. Symmetric 
displacement boundary conditions w=0 were applied 
at the plane of symmetry z=15 mm [w is the nodal 
displacement components along the z-axis (the 
coordinate axes are shown in FIG. 1)]. The models 
were meshed using the 3-D continuous brick finite 
element C3D8R. This element is defined by 8 nodes 
(one at each corner) having three degrees of freedom 
per node (translations in the nodal x, y, and z 
directions). The mesh was condensed towards the 
crack front in order to perform a more accurate strain 
energy release rate analysis. A total of 29 250 elements 
were used in the modeling. Before carrying-out further 
calculations, a mesh sensitivity study was conducted 
with respect to the elements number in order to ensure 
that the mesh was fine enough to give reliable results. 
Typical mesh used in the simulations of the SENB 
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specimen is depicted in FIG. 3. The foam was modeled 
as a linear-elastic material in the finite element 
simulations with modulus of elasticity, E=340 MPa, 
and Poisson's ratio, v=0.300. These elastic properties 
were experimentally measured according to the ASTM 
D 1621-94 methods. The steel impactor was modeled 
as a rigid body because its stiffness is much greater 
than that of the impacted SENB specimen. For this 
purpose, the rigid body modeling option was used. 
All rigid body motions of the impactor were removed 
except the translation in vertical direction. The contact- 
impact algorithm of the ABAQUS/Explict software 
was used to model the interface between the impactor 
and the specimen. In the finite element simulations the 
impactor was initially positioned just above the SENB 
specimen, subjected to gravitational force, and given 
an initial velocity which was estimated for the 
considered drop-height assuming ideal free fall with 
zero initial velocity. 




Time, ms 



FIG. 4 CONTACT FORCE-TIME BEHAVIOUR OF THE SENB 
SPECIMEN WITH NOTCH INCLINATION OF 15° AT IMPACT 
POTENTIAL ENERGY LEVEL OF 9.6 J. COMPARISON BEYWEEN 
THE EXPERIMENTAL DATA (CURVE 1) AND THE FINITE 
ELEMENT COMPUTATIONS (CURVE 2) 

In order to verify the model, the calculated contact 
force-time curve was compared with the measured 
one, as shown in FIG. 4. It can be observed that the 
simulated and the experimental curves are in good 
agreement. Such comparisons were performed for 
various impact potential energy levels at each notch 
inclination used in the experimental testing. Good 
agreement between the computations and the 
measurements allowed us to validate the finite 
element models used to analyze the mechanical 



behavior of the SENB specimens under low-velocity 
impact loading. 

In the present paper the impact fracture toughness 
was characterized in terms of the dynamic critical 

G° 

strain energy release rate, c , at the crack growth 
initiation from the notch tip. 




FIG. 5 SCHEMATIC OF THE CRACK FRONT ZONE WITH 
NODAL NOTATIONS 



The dynamic critical strain energy release rate was 
computed using the virtual crack closure technique. 
The advantage of this technique is the fact that it 
requires only one analysis by the finite element model 
for the actual crack length. The nodal notations in the 
crack front zone are shown in FIG. 5. The virtual crack 
closure technique assumes that the work done to 
closure the crack by one length is equal to strain 
energy released when the crack grows by one element 
length. Therefore, the strain energy release rate mode 

G G G 
components, 7 , 77 , and 777 , at node d in FIG. 5 

can be calculated using the formulae 
1 



2AaAb 
1 

2AaAb 
1 



x u\ u u- u if) 

/ 

Y ld( V le-V lf ) 



(1) 



(2) 



2AaAb t (3) 

where XI, Yl and Zl are the nodal force components, 
and ul, vl, and wl are nodal displacements in the xl, 
yl and zl directions, respectively. The subscripts in 
Eqs. (1) - (3) denote the corresponding nodes in FIG. 5. 

It has to be pointed out that the nodal force 
components XI and Yl in the local coordinate system 
xlylzl, associated with the notch plane (FIG. 5) were 
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computed using the formulae for rotation of 
rectangular coordinates 



X\ = X cos a + Y sin a 



Fj = -X sin a + Y cos a 



(4) 
(5) 



where X and Y are the nodal force components in the 
global coordinate system, a is the angle of notch 
inclination which is defined in FIG. 1. The local 
coordinate system xlylzl is chosen in such way that 
coordinate plane xlzl coincides with the notch plane 
[the zl-axsis is parallel to the crack front, the yl-axis is 
normal to the notch plane (FIG. 4)]. Thus, Zl- and Z - 
nodal force components coincide. The nodal 
displacement components ul and vl were computed 
from u and v using the same formulae for rotation of 
rectangular coordinates (wl and w nodal 
displacement components coincide). 

The components distribution of the strain energy 
release rate mode was obtained applying Eqs. (1) - (3) 
at each node along the crack front. The average values 
of the components of the strain energy release rate 
mode along the crack front were calculated as 



-. w 

Gl=-\G,(z)dz 

rr 

yv 
. w 



G in 



w 



where 



Gr(z) G n (z) 



and 



G m (z) 



are 



(6) 

(7) 

(8) 
the 



components of the strain energy release rate mode 
along the crack front obtained by Eqs. (1) - (3). 

The average value of the total strain energy release 
rate along the crack front was calculated as 



c = g; + G a n + Gj, 



III 



(9) 



The distribution of dynamic strain energy release rate 



Gf G[ 



G D 

mode components, W/ , W// , and 111 , along the 
crack front in the SENB specimen with stationary 
crack was obtained by Eqs. (1) - (3) at each step of the 
transient analysis for the experimentally determined 
critical impact potential energy for the notch 
inclination considered. The average values of the 
dynamic components in strain energy release rate 



^-j Da ^-y Da ^-j Da 

mode, 1 , 11 , and 111 , as well as the average 
value of the total dynamic strain energy release rate 

qDci 

along the crack front, , were calculated by Eqs. (4) 
- (7) at each step of the transient analysis (the integrals 
in Eqs. (6) - (8) were solved numerically). In this 

G^ a 

manner was obtained in a function of the time as 
shown in FIG. 6 for notch inclination of 150 at the 
critical impact potential energy of 16 J. 
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FIG. 6 AVERAGE VALUE OF THE DYNAMIC TOTAL STRAIN 
ENERGY RELEASE RATE ALONG THE CRACK FRONT, G° a , 

PLOTTED AGAINST THE TIME AT THE CRITICAL IMPACT 
POTENTIAL ENERGY OF 14.8 J IN SENB SPECIMEN WITH 
NOTCH INCLINATION OF 15°. THE IMPACT FRACTURE 

TOUGHNESS, G"° = 3970 J/m 2 , IS OBTAINED AT THE 
EXPERIMENTALLY MEASURED FRACTURE TIME, tf =5.5 ms 

The impact fracture toughness (expressed in terms of 
the average value of the total dynamic critical strain 
energy release rate along the crack front at the onset of 

qDo 

crack growth from the notch tip), c , was obtained 
at the experimentally determined fracture time, tf=5.5 
ms, as illustrated in FIG. 6. In this way it was found 
(j Da — 3970 

c J/m2. The distribution of the dynamic 

critical strain energy release rate mode components, 

qD qD qD 

IC , uc , and II1C , along the crack front at tf=5.5 
ms is shown in FIG. 7. Only half of the crack front is 
depicted, because the distribution is symmetrical 
about the crack front centre. The horizontal axis is 
defined such that z/W=0 is at the specimen lateral 
surface. Thus, z/W=0.5 corresponds to the crack front 
centre. It can be observed in FIG. 7 that the 
distribution is non-uniform with strong specimen 
lateral surface effects. Near the specimen lateral 

G, r 



surface, there is a decrease in 



G D 

IC , an increase in ,,c , 
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and appearance of UIC . The mode III component 
arises as a result of different transverse contractions 
which occur in the crack faces. 




0.1 
Lateral surface 



0.4 0.5 z/W 

Mid-plane 



FIG. 7 DISTRIBUTION OF MODE I (CURVE 1), MODE II (CURVE 
2), AND MODE III (CURVE 3) COMPONENTS OF THE DYNAMIC 
CRITICAL STRAIN ENERGY RELEASE RATE AND MODE I 
(CURVE 4), MODE II (CURVE 5), AND MODE III (CURVE 6) 
COMPONENTS OF THE STATIC CRITICAL STRAIN ENERGY 
RELEASE RATE ALONG THE CRACK FRONT IN SENB 
SPECIMEN WITH NOTCH INCLINATION OF 15°. THE 
HORIZONTAL AXIS IS DEFINED SUCH THAT z/W=0 IS AT THE 
SPECIMEN LATERAL SURFACE; THUS z/W=0.5 CORRESPONDS 
TO THE CRACK FRONT CENTRE (W IS THE SPECIMEN WIDTH) 

Such analyses of the components of the strain energy 
release rate mode were carried-out for each notch 
inclination from 150 to 900. Since the components of 
the strain energy release rate mode along the crack 
front was non-uniform, it was decided to characterize 
the fracture behavior by the average values of these 
components which were calculated using Eqs. (6) - (8). 
The corresponding integrals were solved numerically. 



The values of 



GDa i s-% Da 
IC /Cr C 



and the mixed-mode ratios, 

£i Da i Da Da j Da 

lie c ^ an( j mc c r obtained 
were summarized in TABLE 1. The data in TABLE 1 
indicate that the average value of the total dynamic 

qDo 

critical strain energy release rate, c , decreases with 
increasing the notch plane inclination. This is due to 
the transition of the loading conditions in the notch tip 
from mixed-mode I/II at low notch plane inclination 
toward pure mode I at notch plane inclination of 900. 
In other words, when the notch is vertical, the analysis 
generates the pure mode I fracture toughness, while at 
notch tip inclinations less than 900, it generates mixed- 
mode I/II fracture toughness values which are higher 
than the mode I. It is interesting to analyze the 
evolution of the relative amount of mode II 
component of the dynamic critical strain energy 
release rate with increase of notch plane inclination. It 



G ^ a I G ^ a 

can be seen in TABLE 1 that the ratio I!C c 
decreases with increasing the inclination from 0.434 (at 
inclination of 150) to (at inclination of 900). It can be 
summarized that the fracture performance is mode I 
dominated, because at each notch plane inclination the 

G^ a I G^ a 

ratio 10 c is higher than 0.5. TABLE 1 indicates 
also that the relative amount of mode III component is 
very small compared with the mode II and the mode 
III. Therefore, the crack growth in SENB specimens 
containing inclined notch can be regarded as a mixed- 
mode I/II fracture problem. 

The impact fracture toughness was compared with the 
static one at various notch inclinations. The static 
fracture toughness was experimentally characterized 
by the SENB specimens with inclined notch. The 
geometry and dimensions of the SENB specimens 
used in the static testing were the same as ones in the 
low-velocity impact fracture test. The specimens were 
loaded under displacement control at a constant cross- 
head speed of 3 mm/min. The test was stopped after 
about 4 mm crack extension. The critical fracture load 
at the onset of crack growth from the notch tip was 
determined by inspection of the specimen lateral 
surface using a microscope during the testing and by 
observation of the load-displacement response. The 
test data were analyzed by three-dimensional linear- 
elastic finite element models of the SENB specimens. 
The static critical components in the strain energy 
release rate mode (at the onset of static crack growth 
from the notch tip) were calculated by the virtual crack 
closure technique via Eqs. (1) - (3) applied at each 
node along the crack front. At notch inclination of 15o 
the calculations yielded the components of the strain 
energy release rate mode along the crack front that is 
shown in FIG. 7 for comparison with the dynamic 
critical strain energy release rate. The data in FIG. 7 
indicate that the dynamic critical strain energy release 
rate is lower than the static one. The static fracture 

G^ a 

toughness, c , was calculated as the sum of the 
average values of the components of the strain energy 
release rate mode along the crack front using formula 
(9). Such analyses were carried-out for each notch 
inclination from 150 to 900. The obtained results are 
included in TABLE 1. It can be observed in Table 1 
that the impact fracture toughness is between 22.4 % 
(at notch inclination of 150) and 32 % (at notch 
inclination of 600) lower than the static one. This 
finding is attributed to the brittle response of the 
polymer foam to the impact loading. 
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Conclusions 

The impact fracture in rigid foam was investigated by 
SENB specimens with inclined notch. A drop-weight 
rig was used to perform the low-velocity impact 
fracture tests. The impact fracture behavior at various 
notch inclinations was studied in terms of the dynamic 
critical strain energy release rate at the onset of crack 
growth from the notch tip. From the experimentally 
obtained critical impact potential energy and the 
fracture time at various notch inclinations, the 
components of the strain energy release rate mode 
were determined by the virtual crack closure 
technique in conjunction with three-dimensional 
linear-elastic finite element dynamic analysis of the 
SENB specimen using the ABAQUS/Explicit software. 
Based on the obtained results, the following 
conclusions were drawn. 

1. The onset of crack growth from the notch 
tip during the impact testing was signaled 
by a sudden drop in the contact force-time 
response of the SENB specimen. This 
finding is attributed to the specimen 
bending stiffness degradation due to the 
onset of crack growth. 

2. The critical impact potential energy (at the 
onset of crack growth) decreases with 
increase of the notch inclination. This is 
due to the specimen bending stiffness 
decrease with increase of the notch 
inclination. 

3. The developed finite element model can be 
used for analyzing the mechanical 
response of the SENB specimen to impact 
loading, since the computed contact force- 
time curves were in a very good agreement 
with the measured ones at various impact 
potential energy levels and notch 
inclinations. 

4. Non-uniform distribution of the dynamic 
critical strain energy release rate mode 
components along the crack front with 
strong specimen lateral surface effects was 
generated by the analysis. 

5. It was found that the SENB specimen with 
notch inclination varying from 150 to 900 
can be used for experimental investigation 
of impact mixed-mode I/II fracture 
behavior of rigid foams over the range of 

0<G° a c /G° a < 0.434 



6. The impact fracture toughness was 
between 22.4 % (at notch inclination of 150) 
and 32 % (at notch inclination of 600) lower 
than the static one (this finding is 
attributed to the brittle response of the 
polymer foam to the impact loading). Thus, 
the static fracture toughness should not be 
used in the fracture mechanics based 
design of rigid foam structures that are 
subjected to impact loading during their 
lifetime. 
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